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PLAN DU COURS 


OBJECTIF: Conception des algorithmes corrects et efficaces 

PLAN 

□ RAPPELS : NOTATIONS ALGORITHMIQUES 

□ COMPLEXITE 

□ ALGORITHMES ITERATIFS DE TRIS 

□ RECURSIVITE 

□ DIVISER POUR RESOUDRE 

□ PREUVE D’ALGORITHMES 
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Notations algorithmiques 


□ Un type est un ensemble de valeurs sur lesquelles on définit 
des opérations. 

• Types de base : 

y Entier : Opérateurs arithmétiques +, -, *, div, mod 

y Réel : Opérateurs arithmétiques +, -, *, / 

y Booléen : Opérateurs logiques et, ou, non 

y Caractère : constante (lettre imprimable) entre 
apostrophe. 

_ Les opérateurs relationnels permettant de faire des 
comparaisons: <, <, =, >, >, ^ 

Le Résultat de la comparaison est une valeur booléenne. 
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Notations algorithmiques 


□ Une variable possède : 

- un nom 

- une valeur 

- un type 

(la valeur d’une variable peut changer au cours de 
l’exécution) 

Déclaration : <variable> : <type> 


□ Une expression, pour un type, est soit une constante, 
soit une variable, soit constituée à l’aide de constantes, 
de variables, de parenthèses et des opérateurs 
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Notations algorithmiques 
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□ INSTRUCTIONS 

■ Affectation : <variable> := <expression> 


■ Condition : 


si <condition> alors 
action 


fsi 


OU : 


si <condition> alors 

action 1 


sinon 


action 2 


fsi 

( condition est une expression à valeur booléenne; 
action est une instruction ou un bloc d’instructions séparées par ;) 
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Notations algorithmiques 


- Itération 

• boucle Pour 


pour <variable> := <initiale> à <final> faire 
action 

fpour 


(Où initiale et finale sont des expressions de même type que celui de 
la variable contrôlant la boucle, le type peut être entier, caractère 
ou énuméré) 


Remarque: la boucle pour affecte la valeur de initiale à variable et 
compare cette valeur à celle de finale avant d’exécuter action. 
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□ Exemple: calcul de 1 +2 + .. .+n (n entier> 1 
fixé) 

Programme somme_des_n_premiersTermes 
// partie déclaration 

n : entier; s : entier; i : entier; (ou n, i, s : entier) 
Début 

/ / Lecture des données 

Écrire (" n = ? " ); lire(n); 

/ / calcul de la somme 
s := 0; 

pour i := 1 à n faire 
s := s + i; 

fpour; 

/ / affichage du résultat 

écriref'l + 2 + ... + n = ", s); 


fin 


Notations algorithmiques 


en C 


main ( ) { 

// déclaration des variables 
int i, n, s; 

/ / lecture des données 

prinf(" n = ? "); scanf(" %d" , &n); 
// calcul de la somme 
s = 0; 

for(i — 1 ; i <— n; i++) s = s + i; 

/ / affichage 

printff' 1 +2 + ...+ %d = %d \n" , n, s) 

} 
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Notations algorithmiques 


- Itération 

• boucle tantque • boucle répéter 


tantque < condition> faire 
action 

ftantque 


répéter 

action 

jusque < condition> 


Remarques: 

Condition est une expression à valeur booléenne, cette expression 
doit être modifiée dans le corps de la boucle (dans action). 

La boucle pour peut être traduite en boucle tantque (ou en répéter) 
mais l’inverse n’est pas toujours vrai. 
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Programme somme_des_n_premiersTermes 
// partie déclaration 

n : entier; s : entier; i : entier; (ou n, i, s : entier) 
Début 

/ / Lecture des données 

Écrire (" n = ? " ); lire(n); 

// calcul de la somme 
s := 0; 

i:= 1; 

tantque i <= n faire 
s := s + i; 
i := i + 1; 
ftantque; 

// affichage du résultat 

écriref'l + 2 + ... + n = ", s); 


fin 


Notations algorithmiques 


□ en C 
main ( ) { 

// déclaration des variables 
int i, n, s; 

/ / lecture des données 

prinf(" n = ? "); scanf(" %d" , &n); 
// calcul de la somme 
s = 0; 

i — 1; 

while (i <= n) { 

s = s + i; 
i = i + 1; 

} 

/ / affichage 

printff' 1 +2 + ...+ %d = %d \n" , n, s) 

} 
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Notations algorithmiques 


□ Algorithme 

■ Description formelle d’un procédé de calcul qui permet, 
à partir d’un ensemble de données, d’obtenir des 
résultats. 

• Succession finie d’opérations 

DONNEES algorithme, RESULTATS 


Cheminement à suivre : 


Problème 


analyse 


codification 

_ Algorithme _ programme 
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Notations algorithmiques 


D Conception structurée des algorithmes 

Découper l’algorithme en sous algorithmes plus 
simples. Chaque sous algorithme est un algorithme. 

> Fonctions 

Une fonction est un sous algorithme (sous programme) qui, à partir 
de données, retourne un seul type de résultat . 

Une fonction 

• possède un nom 

• communique avec l’extérieur par le biais des paramètres 

• retourne un résultat par l’instruction retourner(expression) 
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Notations algorithmiques 


□ Schéma d’une fonction : 

fonction <nom de la fonction>(Liste des paramètres) : <type du résultat> 
// déclaration des variables locales 

début 

/ / corps de la fonction qui contient l’instruction retourner 

fin 


Les paramètres de la définition d’une fonction (appelés formels) sont : 

• typés 

• séparés par 7 s’il y en a plusieurs 

Les paramètres formels sont utilisés pour ne pas lire les données dans une 
fonction. 

Les variables déclarées dans une fonctions (y compris les paramètres formels) 
sont appelées variables locales. 
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□ Exemple: factoriel n en 
pseudo code: 

Fonction fact(n : entier) : entier 
m, i : entier; 
début 

m := 1 ; 

pour i :—2 à n faire 
m := m * i; 

fpour 

retourner(m); 


fin 


Notations algorithmiques 


□ n! en C : 

unsigned int fact (unsigned int n) { 
unsigned int i, m; 
m = 1 ; 

for (i = 2; i<= n; i++) 
m = m * i; 
return m; 

} 
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Notations algorithmiques 
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o L’appel d’une fonction est utilisé dans une instruction sous la forme : 

<nom de la fonction> (liste des paramètres effectifs) 

> Les paramètres formels et effectifs doivent correspondre en nombre 
et en type 

> Lors d’un appel : 

• les paramètres formels reçoivent les valeurs des paramètres 
effectifs correspondant 

• le corps de la fonction est exécuté jusqu’au premier retourner 
rencontré 

• l’exécution se poursuit (à l’dresse de retour) dans la fonction 
appelante. 

Remarque: le type de retour d’une fonction peut être vide, et dans 
ce cas on écrit retourner)) ou pas d’instruction retourner; 
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□ Exemple: calcul de 1 +2 + .. .+n (n 
entier> 1 fixé) 

Programme somme_des_n_premiersTermes 
// partie déclaration 
n : entier; r : entier; (ou n, r : entier) 
Début 

/ / Lecture des données 

Écrire (" n = ? " ); lire(n); 

/ / appel et utilisation de la fonction 
sommeArith 

r := sommeArith(n); 

/ / affichage du résultat 

écrire)" 1 + 2 + ... + n = ", r); 


fin 


Notations algorithmiques 


□ Fonction pour calculer la somme 

1 +2 + ...+m 



Fonction sommeArith ( m : entier) : entier 
i, s : entier; 
début 
s := 0; 

pour i := 1 à m faire 

S := S + î; 

fpour; 

retourner(s); 

fin 
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Notations algorithmiques 


o Deux types de passages des paramètres : 

• Passage par valeur: la fonction travaille sur une copie du 
paramètre effectif transmis; i.e. la valeur du paramètre effectif 
n’est pas modifiée après l’appel de la fonction. 


• Passage par adresse (ou par référence): 

- l’identificateur du paramètre formel est précédé par le mot ref. 

- la fonction travaille directement sur l’identificateur du paramètre 
effectif; i.e. toute modification sur le formel entraîne la même 
modification sur le paramètre effectif correspondant. 
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Passage par valeur 

Fonction échnger(x: réel, y : réel) : vide 

z : réel; 

début 

Z := x; 

X := y; 

y == z; 

fin 

Fonction appelanteQ 

a, b : réel; 

début 

a := 2; b:= 7; 
échangera, b); 
écrire( a = , a); 
écrire( b = , b); 
fin 


Les résultas affichés par la fonction appelante : 

a = 2 b = 7 


Notations algorithmiques 


Passage par réference 

Fonction echnger(ref x: réel, ref y : réel) : vide 

z : réel; 

début 

Z := x; 

X := y; 

y == z; 

fin 

Fonction appelanteQ 

a, b : réel; 

début 

a := 2; b:= 7; 
échangera, b); 
écrire( a = , a); 
écrire( b = , b); 
fin 

Les résultas affichés par la fonction appelante : 

a = 7 b = 2 
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Notations algorithmiques 



□ Un algorithme se comporte comme une fonction sauf 
que l’on ne s’occupe pas des déclarations des variables 
ni de leurs types. 


> Schéma d’un algorithme : 

Algorithme <nom de l’algorithme> 

Données : // 1 es variables qui sont des données de l’algo. 

<nom de l’algorithme>(Liste des paramètres) Résultat(s): // variable(s) contenant le(s) résultat(s) 


début 

/ /bloc d’instructions Ou 

fin 


début 

/ /bloc d’instructions 

//ne contenant pas retourner 

fin 


On omet la partie déclaration des variables locales en adoptant la 
règle: les variables simples sont en minuscule et les tableaux en 
majuscule. 
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Notations algorithmiques 


□ Algorithme pour calculer n! 

factorîel(n) 

// n > 0 

début 

m := 1 ; i := 1 ; 
tantque i < n faire 
i := i + 1; 
m := m * i; 
ftantque 
retourner(m); 
fin 
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TABLEAUX STATIQUES 


□ TABLEAUX 

■ Un tableau est utilisé pour représenter une suite d’éléments de même type 
( — (t 1 ^2,. . v^ n ))* 

> Déclaration d’un tableau (à un dimension): 

<nom du tableau> : tableau [1 .. max] de <type des éléments> 

(exemple T : tableau [1 ..20] de réel) 

où : - max est une constante entière (positive) 

- le type des éléments est quelconque (de base ou déclaré). 

• La taille (ou longueur) d’un tableau est le nombre d’éléments du tableau 
elle est toujours inférieure ou égale à max. 

• les tableaux, en algorithmique, commencent à l’indice 1 (en C et en java, 
ils commencent à l’indice 0) 
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TABLEAUX STATIQUES 


> Opérations (de base) sur les tableaux: 

❖ Accès à un élément du tableau : 

<nom du tableau> [<indice>] 

(indice est une expression de type entier) 

❖ Parcours : (on utilise un indice et une boucle pour ou tant que) 
(exemple : 

pour i= 1 à n faire //n est le nombre d’éléments du tableau T 

traiter( T[î] ) // traiterQ est une fonction ou algorithme à définir 

fpour; 

❖ Recherche d’un élément dans un tableau 

❖ Insertion d’un élément 

❖ Suppression d’un élément 
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TABLEAUX STATIQUES 
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Exemple : recherche de la position d’un élément dans un tableau de réels. 


fonction localiser(T: tableau[l ..max]de réel, n : entier, val :réel) : entier 
i : entier; 
trouve : booléen; 
début 

i := 1 ; trouve := faux; 
tantque (i < n) et (non trouve) faire 
si (T[i] = val) alors 
trouve := vrai; 
sinon i := i + 1 ; 
fsi; 

ftantque; 

si (trouve) alors retourner(i) 
sinon retourner(O); 
fsi 
fin 
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Algo.ll 


□ insertion d’élément x dans un tableau T à n éléments à la position p. 

Fonction insererfT : tableau[l ..max] de réel, ref n : entier, x : réel, p : entier) : vide 
// 1 <p < n 
i : entier; 
début 
i := n; 

tantque i > p faire 
T[i+l]î= T[î] ; 
i := i - 1 ; 
ftantque 

T[p] î= x; 

n := n + 1 ; 
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Algo.ll 


Suppression d’une valeur d’un tableau T, qui se trouve à la position 
P- 

Fonction supprimer(T : tableau[l .. max] de réel, ref n : entier, p : entier) : vide 
i : entier; 
début 
i := p; 

tantque i < n faire 
T[i] := T[i+ 1 ]; 

i := i + 1 ; 
ftantque; 
n := n - 1 ; 
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Algo.ll 


o Remarques : 

^ Le nom du tableau, dans une liste de paramètres 
formels, est une référence. 

(en C, l’adresse de T[î], notée T + i, est calculée 

adresse(T[i]) = adresse(T[0]) + sizeof(<type des éléments>) * i) 

^ L’insertion (resp. suppression) d’un élément à un 
indice i, nécessite un décalage des (n- î +1 ) 
derniers éléments d’une position à droite (resp. à 
gauche) 
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STRUCTURE (ENREGISTREMENT) 

□ Structures (ou enregistrements) 

■ Le type structure est utilisé pour représenter une suite d’éléments pas 
nécessairement de même type, chaque élément est appelé champs. 

> Déclaration d’un type structure: 

<nom de type structure> = structure 

<variable_champsl > : <type_champsl >; 

<variable_champs2> : <type_champs2>; 

fstructure 

> Déclaration d’une variable de type structure: 

<variable_structure> : <nom de type structure> 
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STRUCTURE (ENREGISTREMENT) 


Par analogie au type struct de C: 

Déclaration : 

struct <nom de la structure> { 

<type_champsl > <variable_champsl > 
<type_champs2> <variable_champs2> 


} 

Déclaration d’une variable de type structure: 
struct <nom de la structure> <variable> 

Ou avec la définition de type : typedef 

typedef struct <nom de la structure> <type_structure>; 

Déclaration d’une variable : 

<type_structure> <variable> 
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STRUCTURE (ENREGISTREMENT) 



□ Exemples : 

Adresse = structure 

numero_rue : entier; 
nom_rue : tableau[l ..20] de caractère; 
code_postal : entier; 
fstructure 

Point = structure 
abscisse : réel; 
ordonnée : réel; 
fstructure 


□ En C : 

struct adresse { 

int numero_rue; 
char[20] nom_rue; 
int code_postal; 

} 

typedef struct adresse Adresse; 

typedef struct point { 

float abscisse, ordonnée; 

}Point; 


Déclaration des variables: 
adr : Adresse; 
p : Point 


/ /Déclaration des variables: 
Adresse adr; 

Point p; 
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STRUCTURE (ENREGISTREMENT) 

> Opération sur les structures : 

❖ Accès à un champs : 

<variable de type structure>. <variable_champs> 
(ex: p.abcisse désigne l’abscisse du poin p) 

❖ Affectation : 

<variable_type_structure> := <variable_type_structure> 

(L’affectation se fait champs par champs) 

Exemple : manipulation des complexes 
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STRUCTURE (ENREGISTREMENT) 


■ Déclaration de type complexe: 

Complexe = structure 

p_reel / p_imag : réel; 

fstructure 

fonction plus(z1 : Complexe, z2 : Complexe) : Complexe 

z : Complexe; 
début 

z.p_reel := zl.p_reel + z2.p_reel; 
z.p_imag := zl.p_imag + z2.p_imag; 
retourner(z); 
fin 


□ Déclaration de type complexe en C: 
typedef struct complexe { 
double p_reel, p_imag ; 

} Complexe; 

Complexe plus(Complexe zl Complexe z2){ 
Complexe z; 

z.p_reel = zl.p_reel + z2.p_reel; 
z.p_imag = zl.p_imag + z2.p_imag; 
retun z; 

} 
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ALGORITHMIQUE II 


NOTION DE COMPLEXITE 



SMI AlgoII 



Complexité 


□ Comment choisir entre différents algorithmes pour 
résoudre un même problème? 

□ Plusieurs critères de choix : 

> Exactitude 

> Simplicité 

> Efficacité (but de ce chapitre) 
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Complexité 


□ L’évaluation de la complexité d’un algorithme se fait 
par l’analyse relative à deux ressources de l’ordinateur: 

o Le temps de calcul 

o L’espace mémoire, utilisé par un programme, pour 
transformer les données du problème en un ensemble 
de résultats. 

L’analyse de la complexité consiste à mesurer ces deux 
grandeurs pour choisir l’algorithme le mieux adapté 
pour résoudre un problème.(le plus rapide, le moins 
gourment en place mémoire) 


On ne s’intéresse, ici, qu’à la complexité temporelle c.à d. qu’au temps de calcul 
(par opposition à la complexité spatiale) SMI A| 9 o11 




Complexité 


□ Le temps d’exécution est difficile à prévoir, il peut 
être affecté par plusieurs facteurs: 

> la machine 

> la traduction (interprétation, compilation) 

> l’environnement (partagé ou non) 

> L’habileté du programmeur 

> Structures de données utilisées 
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Complexité 
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□ Pour pallier à ces problèmes, une notion de 
complexité plus simple, mais efficace, a été définie 
pour un modèle de machine . Elle consiste à 
compter les instructions de base exécutées par 
l’algorithme. Elle est exprimée en fonction de la 
taille du problème à résoudre. 

> Une instruction de base (ou élémentaire) est soit 
une affectation, un test, une addition, une 
multiplication, modulo, ou partie entière. 
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Complexité 


□ La complexité dépend de la taille des données de l’algorithme. 

■ Exemples : 

> Recherche d’une valeur dans un tableau 

— > taille (= nombre d’éléments) du tableau) 

> Produit de deux matrices 
— > dimension des matrices 

> Recherche d’un mot dans un texte 

— > longueur du mot et celle du texte 
On note généralement: 

n la taille de données, T(n) le temps (ou le coût) de l’algorithme. 
T(n) est une fonction de IN dans IR + 
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Complexité 


□ Dans certains cas, la complexité ne dépend pas 
seulement de la taille de la donnée du problème 
mais aussi de la donnée elle-même. 

■^Toutes les données de même taille ne génèrent 
pas nécessairement le même temps d'exécution. 

— » (Ex. la recherche d’une valeur dans un tableau 
dépend de la position de cette valeur dans le 
tableau) 


SMI Algoll 



Complexité 


□ Une donnée particulière d’un algorithme est appelée instance 
du problème. 

□ On distingue trois mesures de complexité: 

1 . Complexité dans le meilleur cas 

T Min (n) — min {T (d) ; d une donnée de taille n } 

2. Complexité dans le pire cas 

T Mx (n) — max {T(cf) ; d une donnée de taille n } 

3. d ans la cas moyen 

T MOY (n) “ 2] p(cf).T(cf) 

d de taille n 

p(d) : probabilité d'avoir la donnée d 
T MIN (") ^ T MOY ( n > ^ T MAx( n > 
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□ Exemple. Complexité de la 
recherche d’un élément x dans un 
tableau A à n valeurs. 


i := i 

tantque (i <n) et (A[i] ^x) faire 
i := i+i; 

ftantque 

si i>n alors retourner faux); 
sinon retoumer(vrai); 

On note par: 

a : le coût d’une affectation 
t : coût d’un test 
d : coût d’une addition 


Complexité 


□ Cas le plus favorable , x est le 

premier élément du tableau: 

T m» = 1 a + 3t 

□ Pire des cas , x n’est pas dans le 
tableau: 

T max( n ) - (n+1 )a +(2n+2)t +nd 

□ En moyenne : 
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Complexité 


■ Complexité en moyenne : 

□ On note par : 

D, (1 <i< n) : ensemble de données (de taille n) où x est 
présent à la i eme position 

D n+1 : ensemble de données où x n’est pas présent 

□ On suppose que la probabilité de présence de x dans une 
donnée est q. De plus, dans le cas où x est présent, on 
suppose que sa probabilité de présence dans l’une des 
positions est de 1 /n 

On a : 

P(Di) = q/n, T(Dj) = i a +(2i+l)t + (i-l)d ; (1 <ï< n) 
p( D n+1 ) = 1 - q , T(D n+1 )=T max = (n+1 )a+(2n+2)t+nd 
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Complexité 



p(Dj) = q/n , T(D.) = i a +(2i+l )t + (i-1 )d ; (1 <i< n) 

P( D n+1 ) — 1 - q, T(D n+1 )— T max (n) = (n+1 )a+(2n+2)t+nd 


T MOY (n) “ 2 p(cf).T(cf) 

d de taille n 

p[d) : probabilité d'avoir la donnée d 

n 

Tmoyin] = Y p(Di)T{Di) + p(Dn + i)T(Pn + 1) 

£=1 


T* 


Tmoyin) = — / (ia + (2i + l)t + (i — 1) d) + (l — q)[(?i 4- 1 )a + (2 n + 2 )t + nd] 

i 


i=l 


vi 4 1 vi — 1 

— ij [ — - — a 4 (n 4 2)t H — d] 4 (l — ç) [(vi 4 l)a+ (2vi 4 2)t 4 nd] 


Cas où q=l , i.e. x est toujours présenttf’moyOrÛ = l / 2 : [(n + 1 )o + (2 n + 4}t + (n - l)rf] 

"SAfii^ot 2tr+dJ. + V 2 (“ + - d) 


On remarque que la complexité de cet algo. est 
de la forme: a n + P , a et P son des constates 
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Complexité asymptotique 
Comportement de T(n) 


□ Pour mesurer la complexité d’un algorithme, il ne s’agit pas de faire un décompte exact du 
nombre d’opérations T(n), mais plutôt de donner un ordre de grandeur de ce nombre pour n 
assez grand. 

□ Notation de Landau 

❖ "grand O" 

T(n) = 0(f(n)) ssi 

3 c > 0 3 n 0 > 0 Vn>n 0 T(n)< c.f(n) 

❖ "grand oméga" 

T (n) = £1 (f(n)) ssi 

3 c > 0 3 n 0 > 0 V n > n 0 T(n) > c.f(n) 

❖ "grand thêta" 

T (n) = 0(f(n)) ssi 

3 c ? > 0 3c 2 >03n 0 >OVn>n 0 c n f (n)< T(n)< c 2 f(n) 


Remarque: les constantes c, c 1# c 2 et n 0 sont indépendantes de n 
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Complexité 



D’une manière générale, f : IR — > IR 

f(x) = O(g)x)) s’il existe un voisinage V de x 0 et une constante k>0 tels 
que |f(x)|< k|g(x)| , (x e V) 


- Si la fonction g ne s’annule pas, il revient au même de dire que le 

. /'CO 

rapport 


£f(0 


est borné pour x G V. 
-Exemple: au voisinage de 0, on a: 
x 2 = O(x), ln(l+x) = O(x) 


- Au voisinage de l’infini (comme pour le cas de la complexité), il 
existe a>0 (V = ]a, + [) et k>0 t.q 

|f(x)|< k|g(x)| , Vx > a 


et on dit que f est dominée asymptotiquement par g. 
(Au voisinage de +°° , on a: x = 0(x 2 ), In x = O(x)) 
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Complexité 



□ Remarques 


1. 

2 . 

3 . 

4 . 


f = O(g) « 
de l’infini. 


le quotient 


est bornée au voisinage 


Il im 

SP — +C-Z 


fi*) I 


a — * f = OÇgJ 


si a — 1, on écrit f ~ g et on a f — 0(g). 

f = 0(g) ne signifie pas que le quotient f(x)/g(x) 
tend vers une limite (1 en particulier). 

Exemple. f(x) = x(2+sin x). On a x < f(x) < 3x 
Vx>0, donc f(x) = 0(x). En revanche, le quotient 
f(x)/x ne tend vers aucune limite lorsque x — > +°o 
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Complexité 


□ Abus de notation: 

On a par définition: 

O(f) = {g: IN— HR/3c>0 3 r> 0 > 0 g (n) < c.f(n) , V n > n 0 } 

£2(f) = {g: IN — HR / 3 c > 0 3n 0 >0 g (n) > c.f(n) , V n > n 0 } 

0 (f) = O(f) n H(f) 

La difficulté, dans la familiarisation avec ces concepts, provient de la convention 
de notation (de Landau) qui veut que l’on écrive : 

g = O(f), ou encore g(n) = 0(f (n)) au lieu de g e O(f) 

De manière analogue, on écrit O(f) = O(g) lorsque 0(f)cz O(g) 

(il en est de même pour les notations £1 ou 0) 
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Complexité 


□ Exemple. Soit la fonction T(n) = V 2 n 2 + 3n 

> T(n) = £2(n) ( n 0 = 1 , c = V 2 ) 

> T(n) = 0(n 2 ) (n 0 = 1, q = V 2 , c 2 = 4) 

> T(n) = 0(n 3 ) (n 0 = 1, c = 4) 

> T (n) = 5 * O(n) 

Supposons que T(n) = O(n) 

3 c >0, 3n 0 >0 : V 2 n 2 +3n < en \/n >n 0 
donc c > V 2 n , contradiction, (la constante c ne peut dépendre de n) 
Remarques 

Si T(n) est un polynôme de degré k alors T(n) = 0(n k ) 

2 . 0(n k ) cz 0(n k+1 ) V K > 0 (idem pour 0) 

0(f(n)) cz 0(f(n)) pour toute fonction f positive 
0(1 ) utilisé pour signifier « en temps constant » 
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Complexité 


□ Remarques pratiques: 

Le cas le plus défavorable est souvent utilisé pour 
analyser un algorithme. 

La notation O donne une borne supérieure de la 
complexité pour toutes les données de même 
taille(suffisamment grande). Elle est utilisée pour évaluer un 
algorithme dans le cas le plus défavorable. 

T(n) < cf(n) signifie que le nombre d’opérations ne peut 
dépasser cf(n) itérations, pour n’importe quelle donnée de 
longueur n. 

Pour évaluer la complexité d’un algorithme, on cherche un 
majorant du nombre d’opérations les plus dominantes. 

Dans les notations asymptotiques, on ignore les constantes. 

L’alqorithme de recherche dans un tableau à n éléments, cité précédemment, est en O(n) 
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Complexité 


□ Ordre de grandeur courant 

0(1) : complexité constante 
0(log(n)) : complexité logarithmique 
O(n) : complexité linéaire 
0(n 2 ) : complexité quadratique 
0(n 3 ) : complexité cubique 
0(2") : complexité exponentielle 
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Complexité 


□ Exemples de temps d’exécution en fonction de la 
taille de la donnée et de la complexité de 
l’algorithme. 

On suppose que l’ordinateur utilisé peut effectuer 
1 O 6 opérations à la seconde (une opération est de 
l’ordre de la |Js) 
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Complexité 






n\T(n) 

log n 

n 

n log n 

n 2 

2 n 

10 

3 \is 

10 |js 

30 |js 

100 \is 

1 000|as 

100 

7 [is 

100 |js 

700 |js 

1/100 s 

10 14 

siècles 

1000 

10 |js 

1 000|is 

1/1 00|is 

1 s 

astrono 

mique 

10000 

13 |js 

1/1 00|is 

1/7 s 

1 ,7 mn 

astrono 

mique 

100000 

1 7 |js 

1/10 s 

2s 
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Complexité 



> Un algorithme est dit polynomial si sa complexité est en 0(n p ). 

^ Un algorithme est dit praticable s’il est polynomial (p < 3). 

Les algorithmes polynomiaux où p > 3 sont considérés comme très lents 
(un algorithme polynomial de l’ordre de n 5 prendrait environ 30 ans pour n=1 000) 
^ Un algorithme est dit exponentiel si sa complexité est supérieure à tout polynôme. 

^ Deux grandes classes de la complexité : 

- IP classe des algorithmes polynomiaux 

- K IP classe des algorithmes « Non déterministe polynomiale » 

On a : 

(0(log n) c O(n) e 0(n log n) e 0(n 2 ) c 0(n 3 ) c 0(2") c 0(e n ) c O(nl) 
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Complexité 


Propriétés 

En utilisant la notation de Landau (pour les fonctions de IN dans IR + ), on a : 
f + O(g) = {f + h/he O(g)} 

(h = f + O(g)oh-fe O(g). 
f O(g) = {f h / h e O(g)} 
f = O(f) 

f = O(g), g = O(h) => f = O(h) 
c O(f) = 0(c f) = O(f) ( c>0) 

O(f) + O(g) = 0(f + g) = 0(max(f, g)) 

O(f) + O(f) = O(f) 

O(f) O(g) = O(fg) 
f = O(g) O g = H (f) 
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Complexité 


□ Calcul de la complexité: règles pratiques 

1. la complexité d’une suite d’instructions est la 
somme des complexités de chacune d’elles. 

2. Les opérations élémentaires telle que l’affectation, 
test, accès à un tableau, opérations logiques et 
arithmétiques, lecture ou écriture d’une variable 
simple ... etc, sont en 0(1 ). 

3. T(si C alors Al sinon A2) = max(T(C),max(T(Al ),T(A2))) 
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Complexité 


4. T(pour i:=e 1 à e 2 faire Ai fpour) = Y t ( ai ) 

i = e 1 

( si Ai ne contient pas de boucle dépendante de î 
et si Ai est de complexité O(m) alors la complexité 
de cette boucle « pour » est 0((e 2 — e n H- 1 )m). ) 

5. La difficulté, pour la boucle tantque, est de 
déterminer le nombre d’itération Nb_iter (ou 
donner une borne supérieure de ce nombre) 

T(tantque C faire A ftantque) = 0(Nb_iter x (T(C) + T(A)) 
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Complexité 


□ Exemples 

î. Calcul de la somme 1 +2+.. .+n 


S: = 0; //0( 1 ) 

Pour i:— 1 à n faire 
S:— S + i; // 0(1 ) 


O(n) 


fpour; 


0(1) + O(n) = O(n) 


T(n) = O(n) 
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Complexité 


2. Calcul de : ^w=v j ‘pour i = 2 , ... ,n 

j=i 

pour i := 1 à n faire 

S := 0; // O(i) 

pour j := 1 à i faire ~ 

S := S + |; // O(i) 

fpour; 

T[i] := s; //O(i) 
fpour; 



T(n) = 0(n 2 ) 
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3. Analyse de l’algo. Suivant 


Donnée n; (n>0) 
Résultat cpt; 

début 
cpt := 1 ; 

tantque n > 2 faire 
n := n div 2; 
cpt := cpt +1 ; 
ftantque 

fin 


Que calcule cet algo? 
Quelle est sa complexité? 


Complexité 


□ cpt = 1 + le nombre d’itérations 

□ Le nombre d’itérations = nombre 
de division de n par 2. 

□ Soit p ce nombre. 

- si n est une puissance de 2, i.e. 
n = 2 P alors p=log 2 (n). 

- p vérifie: 2 P < n < 2 P+1 

p < log 2 (n)< p+1 =>p=E(log 2 (n)) 

□ cpt = 1 + E(log2(n)), cette 
expression de cpt correspond au 
nombre de bits nécessaires pour 
représenter l’entier n. 

□ T(n) = 0(log(n)) 
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TRIS ITERATIFS 


SMI ALGOII 1 

V J 



□ Le problème: 

Etant donnée une suite de n nombres (a lf a 2 ,..., a n ), on 
cherche une permutation (arrangement) des éléments 
de cette suite (a’ lf a’ 2 ,..., a’ n ) telle que a’ 1 <a’ 2 ^...^a’ n . 

♦ A partir de la suite (7,1 ,2,6), un algorithme de tri 
donne comme résultat la suite (1, 2,6,7). 

> On se limite aux nombres entiers rangés dans un 
tableau A à n éléments. 

> Dans le cas où les éléments sont des collections de 
données (enregistrement), on trie le tableau suivant une 
clé (champ de l’enregistrement). 
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□ Tri par sélection: 

Algorithme : 

pour i= 1 à n-1 faire 

-chercher le l eme minimum , A k , de {Â ;+7/ ...,Â n } 
(K est r indice de m/n{A /+ï ,...,A n } dans le tableau A) 

- échanger A t et A k 

L’algorithme fonctionne selon le schéma suivant: 


échange 


A 


trié 





1 i 
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n 
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Analyse du tri par sélection 



□ Algorithme: 

Tri_selection(A,n) 

début 

pour i := 1 à n-1 faire 

/ /Recherche de min{A ; , ..., A n } = A k 
k := i; 

pour /':=/'+ 1 à n faire 

si A[j] < A[k] alors k:=j ; 
fsi; 

fpour 

Il échange de A k et A ; 
temp := A[k ]; 

Mk] := A[i]; 

A[i] := temp ; 

fpour; 

fin 


□ La boucle j détermine le i e minimum; 
elle tourne n-î fois (au maximum) pour 
faire (n-i) tests d’éléments. 

□ Les échanges de A[k] et A[i] 
demandent 3 opérations. 

□ La complexité du corps de la boucle j 
est de la forme 

a(n-i)+b / donc en O(n-i) ( i=l , ..., n-1 ). 

□ La complexité de l’algorithme est de 
l’ordre de: 


n-i 

I 

c = l 


(»-0 


m-l n-l 

= Z fl ~Z ! 

i=l i=l 


n(n — 1 ) 
2 


T(n) = 0(n 2 ) 
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TRI par ISERTION 

□ Algorithme: 

Pour i := 2 à n faire 

- on insère A[i], à sa place, dans le sous-tableau 

A[l..i-1] 

(On cherche le 1 er élément < A[i] parmi {A[i- 1 A[1 ]} en 

décalant d’une position à droite) 

L’algorithme fonctionne selon le schéma suivant: 



i 


i-l i 
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Analyse du tri par insertion 



□ Algorithme : 

Tri_lnsersion( A,n) 
début 

pour i := 2 à n faire 

cle := A[i] ; 
j:= i-1 ; 

Tantque (j>l ) et (cle<A[j]) faire 
A[j+1]:= A[j]; 

i:= i-1; 

ftantque 

A[j+1 ] = cle; 

fpour 

fin 


□ La boucle j tourne, dans le 
pire des cas, (i-1 ) fois 

(i-1 comparaisons et i-1 
décalages) 

□ La complexité du corps de la 
boucle i est de la forme 

a(i-l)+b, pour i=2,...,n. 

□ La complexité de l’algorithme: 


y go 

r=2 


w-l 

ak + & 

k=l 


n(n - 1) 

12- —T 


+ &(n — 1) 


T(n) = 0(n 2 ) 
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Tri à Bulles 



□ On dit qu’on a une inversion s’il existe (i,j) tels que i<j et 


(ai,...,cii, a i+ i,...,cin) ■ =c> (Q^ / ... / Qi_|-i / Qi,... / Q n ) 

□ Un tableau est trié s’il n’a aucune inversion. La complexité du 
tri est proportionnelle au nombre d’inversions qui est de 
l’ordre de c n (nombre de couple(i,j) tels que i<j). 


□ L’algorithme consiste à parcourir le tableau à trier en 
examinant si chaque couple d’éléments consécutifs (a i# a i+1 ) 
est dans le bon ordre ou non, si ce couple n’est pas dans le 
bon ordre on échange ses éléments et ce processus est 
répété tant qu’il reste des inversions à faire. 
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Tri à Bulles 


□ Algorithme!: 

Bulles! (A,n) 

Début 

fini := faux; 

Tant que non fini faire 

i: = 1 ; 

tant que (i < n) et (A[i] <A[i+l] ) faire 
i := i + 1 ; 

ftantque ; 
si i < n alors 



échanger (A[i], A[i+ 1 ]) ; 
fini := faux ; } 

sinon 

fini := vrai ; 

fsi ; 


ftantque ; 



fin 


□ On remarque que les transpositions successives 
font pousser le maximum à la dernière position 
du tableau si on fait un balayage de gauche à 
droite et le minimum à la 1 ere position si on fait 
un balayage de droite à gauche. 

□ Algorithme2: 

Bulles2(A,n) 

Début i := 1; 

tantque i <n-l faire 

pour j := n à i+1 pas - 1 faire 
si A[j- 1 ] > A[j] alors 
éc hanger(A[j], A[j-1 ]); 
fsi; 

fpour; 
ftantque ; 

Fin 

Le nombre de comparaison d’éléments de A < LILli. — 
Le nombre d’échanges < 

T(n) = 0(n 2 ) 
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Exemple 

rw 



Der ech=2 



1 = 2 

Der ech=4 





ALGORITHMIQUE II 


Récurrence et Récursivité 
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Récurrence 


□ Suite récurrente: la définition d’une suite est la donnée 

• d’un terme général défini en fonction du (ou des) terme(s) précédant(s) 

• D’un terme initial qui permet d’initialiser le calcul 

□ Principe de récurrence : 

Soit P un prédicat (ou propriété) sur IN qui peut être soit vrai soit faux (on écrira 
souvent P(n) à la place de P(n) = vrai). 

On suppose que 

■ P(0) vrai 

■ Vn g IN, P(n) => P(n+1) 

Alors , pour tout n e IN, P(n) est vrai. 

Si on considère le prédicat suivant 

P(n) : je sais résoudre le problème pour n 

alors le principe de récurrence nous dit que si je sais résoudre le Pb pour n=0 
et que si je sais exprimer la solution pour n en fonction de la solution pour n+1 
alors je sais résoudre le Pb pour n’importe quel n. 
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□ Examples: 

1 . Puissance 
a 0 = 1 
a n+1 = a a n 

2. Factoriel 
0 ! = 1 

n! = n (n-1 )! , n >1 

3. Suite de Fibonacci 

F 0 = F, = 1 


Ou bien 


n>2 


Récursivité 



i 



a a"' 1 n>0 
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Récursivité 


Un algorithme (ou fonction) est dit récursif s’il est défini en fonction de lui-même. 
Exemples 

Fonction puissancefx : réel, n : entier) : réel 
début 

si n = 0 alors retourner 1 
sinon retourner (x * puissancefx , n-1 )) 
fin 

Factoriel (n) 
début 

si n = 0 alors retourner) 1 ) 
sinon retourner (n*factoriel(n-l )) 
fin 
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■ fact (n) 
début 


fin 


si n — 0 alors retourner(l ) 
sinon retourner (n* fact (n- 1 )) 
fsi 


Récursivité 


□ Le déroulement de l’appel de fact(3): 
3*fact( 

3*2=é \ 2*fact(1 


Fact(3) 

6 



*fact(0) 

* 1=1 


La condition n = 0 est appelée test 
d’arrêt de la récursivité. 


□ Il est impératif de prévoir un test d’arrêt 
dans une fonction récursive, sinon 
l’exécution ne s’arrêtera jamais. 

□ L’appel récursif est traité comme n’importe 
appel de fonction. 
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Récursivité 


□ L’appel d’une fonction (récursive ou non) se fait dans un 
contexte d’exécution propre (pile d’exécution), qui 
contient : 

L’adresse mémoire de l’instruction qui a appelé la fonction 
(adresse de retour) 

Les valeurs des paramètres et des variables locales à la 
fonction. 

L’exécution d’une fonction récursive se fait par des 
appels successifs à la fonction jusqu’à ce que la 
condition d’arrêt soit vérifiée, et à chaque appel, les 
valeurs des paramètres et l’adresse de retour sont mis 
(empilés) dans la pile d’exécution. 
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Récursivité 


L’ordre des instructions par rapport à un appel récursif est important. 

Exemple: 

afficher(n) 

début 


si n > 0 alors 

i « 

afficher(n div 10) 


fsi 


r 


<r~ 


écrirefn mod 10) 


fin 


L’algorithme récursif afficher(n) permet d’afficher les chiffres d’un entier, strictement positif, selon la 
disposition de l’instruction écrie(n mod 1 0): 

- Si l’instruction est placée en I , les chiffres sont affichés dans l’ordre inverse 

- Si elle est placée en T, alors les chiffres seront affichés dans le bon ordre 
Pour n = J 23, on a : 

)—> 3 2 1 

r > i23 


SMI Algoll 



Type de récursivité 


□ Récursivité simple: Une fonction récursive contient un 
seul appel récursif. 

□ Récursivité multiple: une fonction récursive contient 
plus d’un appel récursif (exemple suite de 
Fîbonacci). 

□ Récursivité mutuelle( ou croisée): Consiste à écrire 
des fonctions qui s’appellent l’une l’autre. 

Exemple 
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Récursivité mutuelle 


Pair(n) 

début 

si n = 0 alors 
retourner vrai 

sinon 

retourner (impair(n-l )) 
fsi 
fin 


Impair(n) 

début 

si n = 0 alors 

retourner (faux) 

sinon 

retourner (pair(n-l )) 
fsi 
fin 
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Un peu de Structures de Données 


□ Notion de pile. 

Une pile est une structure pour représenter une suite 
d’éléments avec la contrainte qu’on ne peut ajouter, 
ou enlever, un élément que d’un même côté de la 


ajouter 

som 


pile (dit sommet de la pile). 

□ Exemple pile d’assiettes. 

□ Une pile peut être représentée 

par un tableau et un indice de sommet 


ç 

nlev 
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Notion de Pile 


□ Opérations définies sur les piles: 

initialiser(p : Pile) / /Crée une pile vide. 

sommet(p : Pile) : élément// Renvoie l’élément au 
sommet de la pile p, sous la condition que p soit non 
vide. 

empiler(x : élément, p : Pile) / / ajoute x au sommet 
de la pile p. 

dépiler(p : Pile) / / supprime l’élément au sommet de la 
pile p, sous la condition que p soit non vide . 

pileVide(p : Pile) : booléen / / retourne vrai si p est 
vide. 
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Notion de Pile 


□ Exemple. 

Une expression e est dite bien parenthésée (on se limite au £ (‘ et ‘)’) si : 

î. Le nombre de parenthèses ouvrantes (|e| ( ) est égal au nombre de 
parenthèses fermantes (|e|j) dans e. 

2. Tout préfixe (partie gauche) u de e vérifie: |u| ( - |u| ( > 0. 

Algorithme: on parcourt l’expression e (de gauche à droite). A chaque 

rencontre d’une parenthèse ouvrante on l’empile, et à chaque rencontre 
d’une parenthèse fermente on dépile. 

Si on arrive à la fin de e avec une pile vide, l’expression e est bien 
parenthésée sinon e n’est pas bien parenthésée. 


- l’expression (()())() est bien parenthée. 

- l’expression ())( n’est pas bien parenthésée. 


SMI Algoll 



Transformation du récursif en itératif : 
« Dérécursîvation » 

□ Schéma d’un algorithme récursif: 

algoR(X) 
début 

A 

si C(X) alors 

B; 

algoR(cp(X)) ; 

D; 

sinon 

E; 

fsi; 
fin 


OÙ : 

X : liste de paramètres 

C : condition d’arrêt portant sur X 

A, B, D, E : bloc d’instructions (éventuellement vide) 

Cp(X) : transformation des paramètres 
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Transformation du récursif en itératif 
« Dérécursivatîon » 



□ Algorithme itératif équivalent. 


algoR(X) 


algol(X) 

Début 


p : Pile 

début 


A 

initia liser(p); 

si C(X) alors 

A; 


B; 

tantque C(X) faire 
B • 


algoR(<p(X)) ; 

empiler(X, p); 


D; 

X := Cp (X); 

sinon 


A; 

ftantque; 


• ** 

LU 

E; 

fsi; 


tantque (non pileVide(p)) 

fin 


X := sommet(p); 
dépiler(p); 


ftantque 

fin 
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Dérécusivation 


□ Exemple. 


afficherR(n) 

début 

si n > 0 alors 
a fficher(n div 1 0) 
écrire (n mod 10); 
fsi 
fin 


A = B = E = 0 


afficherl(n) 
p : Pile; 

Début 

initia lïser(p); 
tanque n > 0 faire 
empiler(n / p); 
n := n div 1 0; 
ftantque 

tantque (non pileVide(p)) fai 
n := sommet(p); 
dépiler(p); 
écrire(n mod 1 0); 
ftanque 
fin 
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Transformation du récursif en itératif : 

« Dérécursivation » 

□ Récursivité terminale: 

La récursivité est dite terminale si la dernière instruction 
exécutée est un appel récursive; (Cas où D = 0). Il est 
claire, dans ce cas, d’éliminer la pile dans la version 
itérative. (On dépile pour ne rien faire dans la 2 eme 
boucle). 

□ La récursivité d’une fonction F(X) est aussi dite terminale 
lorsqu’elle se termine par l’instruction retourner(F((p(X))). 
On ajoute, dans ce cas, un paramètre à la liste X pour 
contenir le résultat de retour d’un appel récursive, 
comme le montre l’exemple suivant: 
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Exemple 



fonction FACR(n) ; 

début 

si n= 0 alors retourner (1 ) 
sinon retourner (n* FACR(n-l)) ; 
fin. 

RÉCURSIVE 

fonction FACI(n) ; 
début 


fonction FACR'(n, r) 
début 

si n— 0 alors retourner (r) 
sinon retourner (FACR'(n-7, n*r )) 
fin. 

Ny RÉCURSIVITÉ TERMINALE 

fonction FACI'fn, r) ; 
début 


r: = 1 ; 

tant que n > 0 faire < 

{ r := n* r ; n := n - 1 ; } 

retourner (r) ; 
fin. 

ITÉRATIVE 


tant que n > 0 faire 

{ r : = n* r ; n : = n - 1 ;} 

retourner (r) ; 
fin. 


SMI Algoll 



Complexité des algorithmes récursifs 


complexité des algorithmes récursifs est souvent exprimée 
par une équation de récurrence. 

□ Exemple 1. Complexité de l’algorithme récursif pour 
calculer n! (l’opération dominante est la multiplication) 

Soit T(n) le coût de FACR(n). T(n) vérifie l’équation: 

T(0) = 0 

T(n) = T(n-l) + 1 

la solution de cette équation est : 

T(n) = n = 1 +1 + ... + 1 (n fois) 
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Complexité des algorithmes récursifs 


□ 

□ 


Exemple2. 

Tri par sélection 
sel_rec(T,n) 
début 

si n > 1 alors 

k<— max {i e{l,2,...,n} / T[i] > T[j] , j=l,2,...,n et j ^ i} 

échanger(T[k],T[n]) 

sel_rec(T,n-l ) 


fsi 


/ 


fin 


□ Complexité : T(n) vérifie : 

- T(1 ) = 0 
_ T(n) = T(n-l) + n 

T (n) = n+T(n-l) = n+(n-l )+T(n-2) =...= n+(n-l )+...+2 +C(1) = 


n(n + 1) 
2 


- 1 


T(n) = 0(n 2 ) 
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On ouvre une parenthèse 
Encore un peu de structures de données 

□ Notion d’arbre binaire 


On introduit cette notion pour savoir interpréter les arbres d’appels 
dans le cas d’une récursivité double(où il y a deux appels récursifs). 

- Les arbres sont utilisés pour représenter une suite d’éléments. 

- un arbre est un ensemble de noeuds, chaque noeud représente un 
élément de la suite. 


Définition récursive d’un arbre binaire : 

□ un arbre binaire est: 

Soit vide 
Soit formé : 

■ D’un nœud (appelé racine) 


arbre de racine r 



■ D’un sous-arbre gauche , noté SAG, qui est un arbre binaire) 

■ D’un sous-arbre droit, noté SAD, qui est aussi un arbre binaire 
(les deux sou-arbres gauche et droit sont disjoints). 
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Arbre binaire: terminologie 


Soit A un arbre binaire de racine r. 

□ La racine du SAG (resp. SAD) de A 
est appelé fils gauche (resp. fils 
droit) de r et r est appelé père. 

Un noeud qui n’ a pas de fils est 
appelé feuille. 

un chemin de l’arbre est une suite 
de noeuds n 1/ n 2/ ...n k où n i+1 est un fils 
(gauche ou droit) de ^ , 1 <i<k-l . 

Longueur d’un chemin = nombre de 
noeuds, constituant le chemin, - 1 

Une branche est un chemin de la 
racine à une feuille. 


□ 


□ 


□ 


□ 


□ La hauteur d’un arbre est la longueur 
de la plus longue branche de 
l’arbre. 




-1 si A = 0 


h(A) = 


1 + max(h(SAG(A), h(SAD(A)) 



Branche: 


Branche plus longue 
Hauteur = 4 
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Arbre binaire 


□ Résultat (utile pour la complexité sur les arbres 
binaires de recherche): 


la hauteur h d’un arbre binaire de taille n (n est le nombre de nœuds de l’arbre) 
vérifie: 1 + flog 2 n] < h < n 

- la hauteur est, en moyenne, un 0(log n) et 
un O(n) dans le pire des cas. 


□ Les algorithmes sur les arbres binaires se ramènent 
souvent aux algorithmes de parcours de ces arbres. 
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Algorithmes de parcours 
puis on ferme la parenthèse 


Parcours(a : Arbre) 
début 

si a ^ 0 alors 


parcours(SAG(a)); 


T 




Traiter(racine(a)) 


parcours(SAD(a)); 
T 


<- 


fsi; 


fin 


Selon la disposition de l’instruction 
traiter(racine(a)), on distingue 3 types de 
parcours: 

) : parcours préfixe. 

Y : parcours infixe. 

Y : parcours postfixe. 



On passe 3 fois sur un nœud . 

Parcours préfixe: on traite un nœud lorsqu’on le rencontre pour la 1 er fois. 
10-5-9-50-30-20-40-35- 

Parcours infixe: " " " “ la 2è fois. 

5-9-1 0-20-30-35-40-50-60 

Parcours posfixe: " " " le quitte pour la dernière fois. 

9-5-20-35-40-30-60-50- 
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Récursivité double 

& 

Arbre des appels récursifs 


□ Exemplel : calcul de la hauteur 
d’un arbre binaire. 

h(a : arbre) 

début 

si a = 0 alors retourner -1 
sinon 

hl := h(SAG(a)); 
h2 : h(SAD(a)) ; 
retourner) 1 +max(h1 ,h2)); 
fsi; 
fin 


□ Un arbre est donné par sa racine 

□ appel de h(l ): 
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Récursivité double 

& 

Arbre des appels récursifs 


Exemple2 : tours de Hanoi (Occupation des moines de Hanoï) 

Le jeu consiste à faire passer les 
disques de la tour Tl à la tours 
T2, en ne déplaçant qu’un seul 
disque à la fois, et en utilisant la 
tour intermédiaire T3 de telle 
sorte qu’à aucun moment un 
disque ne soit empilé sur un 
disque de plus petite dimension. 

La solution semble difficile, et pourtant une solution récursive existe. 

Soit n le nombre de disques à déplacer. Si n=l la solution est triviale. 

Si on sait transférer n-1 disques alors on sait en transférer n. 

Il suffit de transférer les n-1 disques supérieurs de la tours Tl vers la tours T3, 
de déplacer le disque le plus grand de Tl vers T2, puis de transférer les n-1 disques 
de T3 vers T2. Ceci se traduit par l’algorithme récursif suivant: 
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Tours de Hanoï 



H(n,Tl ,T2,T3) 
début 

si n = 1 alors écrire(Tl, >’,T2) 
sinon 


H(n- 1,11,13,12); 
écrire(Tl /— »’,T2) 
H(n-1,T3,T2,T1); 
fsi 
fin 
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H(n,Tl ,T2,T3) 
début 

si n = 1 alors écrire(Tl, >’,T2) 
sinon 


H(n-1 ,T1 ,T3,T2); 
écrire(Tl >’,T2) 


H(n-1,T3,T2,T1); 

Arbre des appels de H(3,a,b,c) 


H (3, a, b, c) 



H (2, a 

, c, b) 

a -> b 

H(2,c 

:,b,a) 







-> b 


H(1 ,a,b,c) 

a — » c 

H(1 ' b ' c ' a ) H(1 ,c,a,b) 

c 

H(1 ,a # b,c) 


a — » b 


b — » c 


c — > a 


a -» b 
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Tours de Hanoï : Complexité 



Soit T(n) le temps pour déplacer les n disques. T(n) vérifie l’équation 


“T(l)=l 

T(n) = 2 T(n-l) +1 


On a : 

T(2) = 2+1 

T(3) = 2(2 +1) + 1 = 2 2 + 2 + 1 
On montre, par récurrence, que 

T(n) = 1 + 2 + ... + 2 n_1 = 2 n - 1 

Sachant que T(1 0) = 2 10 —1= 1 023 et une année = 0.3 xl O 8 secondes, il 
faudrait, pour les moines, 10 10 siècles pour pouvoir déplacer 64 disques! 
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Méthode « Diviser pour Résoudre » 


E. CHABBAR 
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DVR 


□ Les algorithmes sont regroupés en familles selon 
certains concepts t. q. Division pour Résoudre, 
Glouton, Programmation dynamique. 

□ - L’aspect DVR consiste à diviser le problème initial 
(de taille n) en sous-problèmes similaires de tailles 
plus petites (en général de taille n/2, n/3,...). 

- Ces sous-problèmes sont résolus récursivement. 

- On peut combiner ces solutions récursives pour 
avoir la solution du problème initial. 
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DVR 



Calcul de x n (n > 1 ) 

Le nombre de multiplications par la méthode classique ( x n = x x n_1 ) est de l’ordre 
de n. 

La méthode DVR exploite la définition suivante: 

1 si n=l 

x n = x p . x p si n= 2p 

x . x p . x p si n = 2p+l 

- pour calculer x n on fait appel (récursif) au calcul x p ou p = E(n/2). (algo. Slide 
suivant) 

- si T(n) est le nombre de multiplications pour calculer x n alors T(n) = T(n/2) + 1 si n 
est pair ou T(n) = T(n/2) + 2 si n est impair. On a, dans tous les cas, T(n) = T(n/2) + 
0 ( 1 ). 

Comme il y a log 2 n divisions euclidiennes successives (avant d’avoir 1 comme 
quotient) , donc T(n) = 0(log(n)). (pour s’en convaincre prendre n une puissance de 
2 ) 

(T(n)=0(l )+0(l )+...+0(l ) log 2 n fois) 
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DVR 

Exemples 


i. Calcul de x n (n > 1 ) 

puissBin(x,n) 

Début 

si n = 1 alors retourner x 
sinon 

si n mod 2 = 0 alors 

y := puissBin(x / n/2) 
retourner (y*y) 

sinon 

y := puissBin(x / n/2) 
retourner (x*y*y) 
fsi; 
fsi; 

Fin. 
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DVR 

Exemples 


2. RECHERCHE DICHOTOMIQUE 

□ La recherche séquentielle d’un élément x dans un 
tableau A à n éléments est de l’ordre de n. 

□ La recherche dichotomique exige que le tableau 
soit trié. Elle consiste à: 

- comparer x à l’élément du milieu 

- si c’est différent, x peut se trouver soit dans la 
moitié gauche soit dans la moitié droite du 
tableau A, selon que x < A[milieu] ou x> A[milieu]. 
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DVR 

Exemples 

2. Recherche dichotomique 

rechDichofAjinfjSupjX) 

début 

si inf < sup alors 

m := (inf + sup)/2 
si x = A[m] retourner m; 

si x < A[m] alors retourner rechDichofA^n^m-l ,x); 
sinon retourner rechDicho(A,m+l ,sup,x); 
fsi; 
fsi; 

sinon retourner 0; 
fsi; 
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DVR 

Exemples 


2. Recherche dichotomique 

□ Complexité de la recherche dichotomique 

Dans le pire des cas (i.e. x n’est pas dans le tableau) 
le nombre de comparaisons T(n), pour chercher x dans 
un tableau A[1 ..n] à n éléments, vérifie: 

T(l)= 1 

T(n) = T(n/2) + 1 , n>l 

- la solution de cette équation dépend du nombre 
d’itérations (nombre de divisions par 2) 
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DVR 

Exemples 


Itération 

0 

1 

2 

3 


Nbre d’élts du ss tableau 
n 

n/2 

n/4 

n/8 


p n/2 p 

L’algorithme utilise p itérations(et aussi p comparaisons) et s’arrête 
lorsque n/2 p =1, c.à d. p = log 2 n. 

par conséquent, la recherche dichotomique est en 0(log 2 n), i.e. 

T(n) = 0(log 2 n) 
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3. Tri par fusion 

■ Fusion de deux tableaux triés. 

Etant donnés deux tableaux triés Tl [1 ..ni ] et 
T2[l ..n2]. La fusion consiste à construire un tableau 
T[1 ..ni +n2] contenant tous les éléments de Tl et 
T2 dans l’ordre croissant. 


Tl T2 




Fusion 



Fusion(Tl ,nl ,T2,n2) 

//T est un tableau qui contient le résultat de la fusion 
il : — 1 ; i 2 : — 1; k:=l; 
tantque (il <nl ) et (i2<n2) faire 
si Tl [il ] < T2[i2] alors 
T[k] Tl [il ]; 
k: = k+l ; il il +1 ; 

sinon 

T[k] := T2[Ï2]; 
k: — k+ 1 ; i2 := i2 + 1; 
fsi; 

ftantque; 

// on recopie les élts restants dans l’un des //tableaux Ti dans le tableau T 

si il < ni alors copier(Tl ,il ,nl ,T,k) 
sinon copier(T2,i2,n2,T / k) 
fsi; 

retourner (T); 


copierfA^f^j) 

début 

pour i:— d à f faire 

m := A[i] 

i i+1; 

fpour; 

retourner(B); 

fin 
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la complexité de la fusion est en 0(n1 +n2) 




TriFusion 


□ Le tri par fusion est un exemple typique de la 
stratégie diviser pour résoudre, à savoir: 

î. Diviser le tableau T[1 ..n] en deux sous-tableaux 
T[l..E(n/2)] et T[E(n/2)+l ..n] 

2. Trier (récursivement) les deux sous-tableaux (2 
appels récursifs à la même fonction TriFusion). 

3. Fusionner les deux sous-tableaux; 

Ceci est schématisé par l’exemple suivant: 
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TriFusion 


Diviser 


Résoudre récursivement 


Combiner 
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TriFusion 



TriFusion(T, n) 

// T 1 # T2 : des tableaux (locaux) de longueurs variables à chaque appel 

début 


si n > 1 alors 

copier(T, 1 , n/2, Tl ,1 ); 
copier(T, n/2 + 1 , n, T2, 1 ); 

Tl := TriFusion(Tl , n/2); 

T2 := TriFusion(T2, n — n/2); 

T := Fusion(Tl, n/2, T2, n - n/2); 


fsi; 

retourner(T); 



fin 


Retourner(Fusion(TriFusion(T 1 ,n/2),n/ 2,TriFusion(T2 / n-n /2) / n-n / 2) 
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TriFusion: exemple d’exécution 

Tableau à trier: 


3 

9 

1 

6 

5 

2 

4 

8 



TF(1 ..4) 


3 

9 

1 

6 


TF(1 ..8) 



1 ..4,5. .8) 


8 



TF(l..l) TF(2..2) F( 1 -.2) 


TFÜ 

~2) 


TF(3 

■Æ 

FM 

..2,3 

-.4) 

3 

9 


1 

6 

1 

3 

6 

9 


TF(3..3) TF(4..4) F(3..4) 


1 


TF(5..6) 


TF(7..8) 


8 


TF (5..5) TF(6..6) F(5..6) TF(7..7) TF(8..8) 



i \ . 

2 

4 

5 

8 

) F(7..8) 


4 

8 



Recopie de ss-tableau 



Tableau (vert)retourné 
à l’appel récursif 


Fusion des 2 ss-tableaux(rouges) 
résultats des appels récursifs 
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DVR: Compléxité 

□ La complexité T(n) pour trier un tableau de n 
éléments par l’algorithme TriFusion vérifie: 

/ T(1 ) = 0 

T(n) = 2 T(n/2) + O(n) , n > 1 


(Il y a 2 appels récursifs, chacun porte sur la moitié du 
tableau. O(n) pour recopier les 2 ss-tableaus en 
2xO(n/2) + leur fusion en O(n))). 
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Equation de récurrence des alg.DVR 


□ La récurrence, utilisée par les algorithmes type DVR, est 
souvent de la forme: 

0(1 ) pour n=l 



a T(n/b) + 0(n d ) , n > 1 


a : est le nombre de divisions du problème initial en sous-problème 
(nombre d’appels récursifs) 


n/b : la taille de chaque sous-problème (b>2) 


0(n d ) : le temps nécessaire pour décomposer le problème en sous- 
problème + le temps pour combiner les solutions des ss-problèmes 
pour avoir la solution du problème initial. 
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DVR: Résultat de Complexité 


□ Théorème: 

Soit T : IN IR+ une fonction croissante à partir d’un certain rang, telle 
existe des entiers n 0 >l, b>2 et des réels d>0, a>0 pour lesquels 


T(n o y - te 
T(n) = aT + en* 


Alors on a : 


T(n) = 


O (n d ) 


0(71 d logfTl ) 

O (n logba ) 


si a 


71 

n > îIq, — puissance de b 

7Lq 


< b d 


SI 


i a = b d 


si a > b d 
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DVR: Résultat de Complexité 


soif T(n\ = OT ('^ j + cn d , 

On montre par récurrence sur p (p>l ), que : ^00 = + cnd ^ (^V 


y-i 


i=0 


prenons — = b? ,et donc p = OQog b n).Larelationprédente devient t 
n Q 

L p - i 

F(jï) = -I- c?i d ^ (^) J ' = /(îî) + cn d ^(n) , ot : 

/=o 

f(n) = ka ? =0(n lo m a ) et g{n) = ZSS) 7 


dUi/cutfruffs a p = e ?lna = e loë& »c lnc; = e“ In b = (— ) Io ® a = 0(n l0Ê:£?a ) 

a ^ f"P— i 

g{n)estune suite géométrique déraison r = — , , par suite g(n) = (r 1) 


Irr 


1 1 j t 


n 


b ü 


r— 1 


premier cas: sir= = 1, alors g{n} = p= O(log & n) et par conséquent = o[n d logiîl} si 0 = b d 

deuxième coj :g{n) = — < sir = — r < 1 et g (V) = 0(1) dans ce cas , p&r conséquent rOl) = sï fî < Ô £ 

1 — r 1 — r b a 


^ 

troisième cas-.g(n) = ■< r p sir = . d >\etg(yi) = 0{îi -£i /^aïre- 


a 


b d 


d'où T[n) = o( n loEîa } sï a > b' 


a? 

r'P — 

b* d 

SMI-ALGOII 


"Jt! 

= 0(w lo ^ E )(— )~ d = 0(n lo ^ a )0(n- d ) 
n a 




DVR: Résultat de Complexité 


D’une manière générale. Si 

T(1 ) = 0(1 ) 

T(n) = a T(^ ) + 0( n rf ) ( n > 1 ,a> 0 , b> 1 , d > 0) 

Alors 


0(n d ) 


si a < b d 



0(n d logera ) 




si a = b d 
si a > b d 
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DVR 


□ Applications 

î. Pour la recherche dichotomique et le calcul de la 
puissance on a : T(n) = T(n/2) + 0(1 ) 

donc T(n) — 0(log n) (a= 1 , b=2, d=0) 

2. Pour le TriFusion on a : T(n) = 2 T(n/2) + O(n) 

Alors T(n) = 0(n log 2 n) (a=2, b=2, d= 1 ) 

s. T(l) = 1 

T(n) = 2t(n/2) + 0(n 2 ) 

a pour solution T(n) — 0(n 2 ) (a=b=d=2) 
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DVR : tri rapide (quickSort) 


□ Soit à trier le tableau T[g..d] (au départ g— 1 et d— n) 

□ Le principe de l’algorithme réside dans une procédure, 
appelée partition, qui réorganise les éléments de T 
autour d’un pivot (élément du tableau choisi au hasard) 
de sorte que : 

1 ) Il existe un indice p (g< p <d) tel que p est la 
position définitive du pivot (T[p] — pivot). 

2) tous les éléments T[g], ..., T[p-1] sont inférieurs ou 
égaux à T[p]. 

3) tous les éléments T[p+1], ...T[d] sont supérieurs ou 
égaux à T[p]. 
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DVR : tri rapide 



□ Le travail de la fonction partition consiste à: 

o Choisir un élément du tableau comme pivot(par exemple le premier 

T[g]) 

Parcourir le tableau depuis la gauche(de gauche à droite) jusqu’à 
rencontrer un élément > T[g] 

Parcourir le tableau depuis la droite (de droite à gauche) jusqu’à 
rencontrer un élément < T[g] 

o Echanger ces deux éléments dans le tableau 

o Continuer ce processus jusqu’à ce que les deux indices (de gauche et 
de droite) se croisent. 


o 


Finalement, é changer le piv ot T[g] et l’élément indiqu é par l’in dice 
de droite. 


< pivot 


Non encore analysé 


> pivot 
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► i 


I 


g 


d 




DVR : Tri Rapide 


/ /T[n+1 ] = + oc 
Partition(T, g, d) 
début 

pivot := T[g]; 

i g; i := d+ 1 ; 

tantque i < j faire 

i := i +1; 

tantque T[i] < pivot faire i:=i+l; 
ftantque 

i := ï -1; 

tantque T[j] > pivot faire j:=j-l; 
ftantque 

si i < j alors échanger(T, i, j); fsi 
ftantque 
échanger(T, g, j); 
retourner(j) 
fin 


□ Exemple 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (1 0) (1 1 ) (1 2) (1 3) 
3141592653 5 8 9 ij 

3141592653 5 8 9 3 10 

3131592654 5 8 9 3 10 

3131592654 5 8 9 57 

3131295654 5 8 9 57 

3131295654 5 8 9 65 

2131395654 5 8 9 5 


1 3 1 


9 5 6 5 4 


8 
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DVR : Tri Rapide 

□ La fonction partition, appliquée à un tableau T, 
produit trois sous-tableaux: 

o Un sous-tableau réduit à un seul élément T[p] qui 
garde sa place définitive dans le tri de T, et 

o Deux sous-tableaux T[g .. p-1], T[p+1 .. d]. 

□ Pour trier T, il suffit d’appliquer récursivement le 
même algorithme sur les deux sous-tableaux. 
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DVR : Tri Rapide 

quickSort(T, înf, sup) 
début 

si inf < sup alors 

p := partition(T, inf, sup); 
quickSort(T, înf, p-1 ); 
quickSort(T, p+1, sup); 
fsi 
fin 
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Complexité du Tri rapide 


La complexité de la fonction partition, appliquée à 
T[1 ..n], est en O(n). 

- Cas le plus défavorable : 

Cas où le pivot sort, à chaque fois, en premier 
(ou en dernier) élément (T: tableau trié). 

La partition coupe le tableau en un morceau de un 
élément et un morceau de n-1 éléments, dans ce 
cas on a : 

C(n) = C(n-l) + O(n) 

(O(n) est le coût de la partition) 

On en déduit que C(n) est en 0(n 2 ) 
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Complexité du Tri rapide 

- Cas le plus favorable : 

cas où le pivot est l’élément médiane de T. 

La partition coupe T en deux morceaux de 
taille n/2 

C(n) = 2 C(n/2) + O(n) 

ce qui donne: C(n) = 0(n log n) 

La complexité moyenne est aussi de l’ordre de n log n 
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Complexité du Tri rapide 



- Complexité moyenne du tri rapide 

La formule de récurrence donnant le nombre de 
comparaisons effectuées par le tri rapide pour une 
permutation aléatoire de n éléments vérifie 

C 0 = C, = 0 et 


i 


n-l 


C n = "-1 + - I(C, + C 

i=0 


pour n >2 

Le terme générique C n peut s’écrire : 


o n-l 

C n = " - 1 + TP I Ci 

i=0 
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Complexité moyenne du tri rapide 


n-1 


• C n n 1 + n î^D C ' 


n-2 


C n = n — 1 + — C , + ^ 

n n i=0 


IC, 


C , + 

I 


n-2 


C — n — 1 + v “'n-l 

n n 


n-1 o 

(n-2+ 2 Z Q) 

n-1 i=D 


n 


( n-1) (n-2) 
n 


C n = 


c = 


n 


C n + 

n-1 


n-1 


n 


C i + 

n-1 


2n-2 


n 


n+1 


n 


C , + 

n-1 


2(n-l ) 


n 


En posant D n = 


n+1 


On aura la récurrence : 


D — D n + 

n n-1 


n+1 n(n+1) 


En négligeant le dernier terme de D n on a: D n — log n 

Du fat que : 1 + ~T~~ + ••• + 1 

n 


— - c " est en 0(n log n) 




Preuve de programmes 


Notions de Logique 


E. CHABBAR 
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Logique 


La logique a joué un rôle décisif dans le 
développement de l’informatique, notamment en 
informatique théorique: 

- Définition d’un modèle théorique des premiers 
ordinateurs (Machine de Turing) 

- Calcul booléen pour la conception et l’étude des 
circuits. 

- La récursivité pour définir la calculabilité. 

- la décidabilité et la complexité pour étudier la limite 
de la machine. 

- La programmation fonctionnelle. 
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Logique 


□ Actuellement l’informatique a envahi tous les 
domaines de la vie. Le problème le plus important 
qui se pose dans la conception d’une application 
informatique est de prouver qu’elle est exempte 
d’erreurs et qu’elle résout le problème pour 
laquelle elle a été conçue. Pour cela, on définit une 
tache par une formule de logique et on montre 
qu’elle est satisfaite par un modèle de cette 
application. Ce type de preuve est appelé « 
Vérification Formelle ». 
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Logique 


□ Une logique, par définition, est un ensemble de formules. 

□ Une formule est construite, sur un alphabet, suivant certaines 
règles (Syntaxe). 

□ La sémantique d’une formule est une valuation (ou 
interprétation dans un modèle) qui détermine la valeur de 
vérité de la formule. 

□ Il y a plusieurs types de logiques: 

- Logique des propositions (d’ordre 0) 

- Logique du premier ordre (les prédicats en font partie) 

- logique du second ordre et logique d’ordre supérieur. 
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Logique des propositions 


□ - On note P = {p, q, ...} l’ensemble des propositions 
atomiques. Chaque proposition atomique est une 
variable qui ne peut prendre que « vrai » ou « faux ». 

{ 1 / a, v t •£=>} l’ensemble des connecteurs (ou 
opérateurs) logique plus les parenthèses. 

□ - On note F(P) l’ensemble des formules déduit de P. 

□ F(P) est le plus petit ensemble qui contient P et qui est 
stable par les connecteurs. 

En d’autres termes: 

Une formule est une suite de symboles de Pu {-i , a, 
v, =>, <=>, (, )} construite selon les règles suivantes: 
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Logique des propositions 


- Toute proposition atomique est une formule. 

- si f est une formule alors — if est aussi une formul 

- si fl et f2 sont des formules alors : 

- fl a f2 est une formule 

- fl v f2 (( « « 

- fl => f2 

- fl <=> f2 « (( (( 

Exemple: (— ip => q) v (p a q) est une formule. 
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Logique des propositions 


□ La sémantique des opérateurs logiques est donnée 
par des tables de vérité. 
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Prédicats 


□ Les prédicats sont construits avec : 

- les constantes (0,1,2,....) 

- les variables (x, y, ....) 

- les fonctions (f, g, +, *, ...) 

- des connecteurs logiques (— i , A, v, =>, 

- des parenthèses 

- des quantificateurs (V, 3) 

□ Un prédicat atomique est un prédicat qui ne contient ni 
connecteur ni quantificateur. 

- Exemples: x < y ; pair(2 x) ; 
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Prédicats 


Règles de formations des formules de prédicats 

- Tout prédicat atomique est une formule 

- si fl et f2 sont des formules alors 

— i f 1 , f 1 v f2, fl Af2, fl => f2, f 1 <=> f2 
sont des formules. 

- si f est une formule alors 

V x f est une formule 
3 x f est une formule 
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Prédicats 


□ Exemples: 

1 ) Vx pair(x + x) 

2) 3x premier(x) a x < succ(succ(0)) 

3) Vx (oiseau(x) => vole(x)) (tous les oiseaux volent) 

4) 3x (oiseau(x) a — i vole(x)) (4 = — i 3) 

5) Vi Vj (1 <i <j<n — ^ T[î] ^ T [ j] ) (spécification d’un tableau 

□ Remarques 

-i (Vx P) = 3x — iP 
-i (3x P) = Vx — iP 
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Déduction logique 


□ Démonstrations logiques 

□ Un séquent est un couple de la forme (ÿ ,f), où f 
est une formule et $ un ensemble fini de 
formules . L’ensemble $ est l’ensemble 

des prémisses du séquent, la formule f sa 

conclusion. 
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• Séquents prouvables 

• Un séquent ,f) est prouvable , ce que l’on notera 
ÿ \- f, s’il peut être construit en utilisant un nombre fini 


de fois les 6 règles suivantes 


1 . 

2 . 

3 . 

4 . 

5 . 

6 . 


si fe^, alors ÿ |-f 

si g et ÿ |-f, alors ÿ, g |- f 

si $ \- (f => f’) et ÿ f, alors f \- f ’ 

si ^ , f |-f’, alors 

^ |-fssi^ h— 1 f 

si ^ ,f |- f’ ÿ ,f |- -if alors $ |— 'f 


(Hypothèse) 


(Modus Ponens) 
(Synthèse) 


(Raisonnement par l’absurc 
(Contradiction) 
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Démonstrations logiques 

Une démonstration d’un séquent prouvable ÿ |- f est 
une suite finie de séquents prouvables f\ |- fi, 

i = 1 ,...,n, telle que : 
j n = j et f n = f 

chaque séquent de la suite est obtenu à partir des 
séquents qui le précède en appliquant l’une des 6 règle 

Remarque : le premier séquent de la suite est 
nécessairement obtenu par utilisation d’une hypothèse 
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• Démonstrations logiques - Exemple 

• Tous les hommes sont mortels, et 

• les Grecs sont des hommes, donc les 
Grecs sont mortels 

h = « être un homme » 

m = « être mortel » 

g = « être Grec >> 

(h m), (g => h) |- (g => m) 
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• Démonstrations Logiques - Exemple 

• Tous les hommes sont mortels, et 

• les Grecs sont des hommes, donc les Grecs 
sont mortels 


1 . 

(h => m), (g => h), g 

-g 

hypothès 

2. 

(h => m), (g => h), g 

- (g => h) 

hypothès 

3. 

(h => m), (g => h), g 

-h 

MP (1&2) 

4. 

(h => m), (g => h), g 

- (h => m) 

hypothès 

5. 

(h => m), (g => h), g 

- m 

MP (3&4) 

6. 

(h => m), (g => h) ^ (g => m) 

synthèse 
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Preuve de programmes 




Logique de Hoare 
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Correction de programmes / 

Spécification 


□ un programme est correct s’il effectue la 
tâche qui lui est confiée dans tous les cas 
permis possibles 


□ nécessité de disposer d’un langage de 
spécification permettant de décrire 
formellement la tâche confiée à un 
programme. 
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Correction de programmes / 

Spécification 

□ description des propriétés que doit 
satisfaire un programme pour répondre au 
problème posé 

- relation entre les entrées et les sorties 
du programme 


Entrée E — 
Précondition : P(E) 


Programme 


— * Sortie S 

Postcondition : Q(E, S) 
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Correction de programmes / 

Spécification 


□ Exemple 

une spécification pour le problème du 
calcul de la racine carrée entière par 
défaut : 

Données n : entier ; 

Résultat r : entier ; 


Pré condition : n > 0 

Algorithme: r := 0; tantque (n >(r +1) 2 ) faire r := r + 1 ftantque 
Post condition : (r 2 < n) a (n < (r+1) 2 ) 


SMI-ALGOII 







Test vs Vérification 



• Les couples (n, r) de l’ensemble suivant {(0, 0), (1, 1), (2, ), (3, 1), (4, 2), ...} 
sont des tests qui réussissent 

• Hoare propose une exécution de l’algorithme sur une valeur 
symbolique qui est un ensemble de valeurs défini par une expression 
logique (ou prédicat). 

• L’ensemble {(0, 0), (1, 1), (2, 1), (3, 1), (4, 2), ...} est 

défini par {(n,r) / n>0 a r >0 a r 2 < n a (n < (r+1) 2 } 

• On note {n>0 a r >0 a r 2 < n a (n < (r+1) 2 } la valeur symbolique des 
variables n et r. 

• L’exécution d’un algorithme A sur une valeur symbolique de données 
définies par l’expression { p} qui donne une valeur symbolique résultat 
définie par l’expression { q} est notée {p} A {q} . {p} A {q} est appelé triplet 
de Hoare (p : précondition, q : postcondition). 
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Système formel 


• Interprétation d’un triplet de Hoare 

{p} A {q} signifie : 

Si la propriété p est vraie avant l’exécution de A ET si l’exécution de A 
se termine, ALORS la propriété q est vraie après l’exécution de A. 
(Correction partielle) 

• Un système formel est un triplet <L, Ax, R> où : 

- L est un langage définissant un ensemble de formules, 

- Ax est un sous-ensemble de L ; chaque formule de Ax est appelée axiome, 

R est un ensemble de règles de déduction de formules à partir d’autres formules 
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Langage algorithmique 

Instructions : 

- Affectation (:= symbole d’affectation et = pour la comparaison) 

- Contrôle : 

Si <cond> alors instruction fsi 
ou Si <cond> alors instruction 
sinon instruction fsi 

- Boucle : 

Tantque <cond> faire instruction ftantque 
Restriction : 

- Pas de fonction (ou procédure) 

- ni de variable pointeur 

- pas de désignation de la forme t[i] où i fait référence à un autre 
tableau. 
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Logique de Hoare 

• Définition de la logique de Hoare : La logique de 
HOARE est un triplet <L, Ax, R> avec : 

- L est l’ensemble des formules { p} A { q} où p 
et q sont des prédicats et A est un fragment 
d’algorithme(ou de programme) 

- Ax est l’ensemble des axiomes de la forme 
suivante : { p[ x / e] } x := e { p(x)} 

(p(x) est obtenu de p[x/e] en substituant toute 
occurrence de e par x) 

- R est l’ensemble de règles de déduction 
suivant : 
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Logique de hoare 


(Séq) si {p} Al {r} ; {r} A2 {q} alors 

{p} Al ; A2 {q} 

(cons g ) si p p’, {p’} A {q} alors 

{P} A {q} 

(cons d ) si {p} A {p’}, p’ q alors 

{p} A {q} 

(condj si {p A c} a {q} , (p A — ic) =^> q alors 

{p} si c alors A fsi {q} 

(cond 2 ) si {p A c } Al {q},{p A — ic} A2 {q} alors 

{p} si c alors Al sinon A2 fsi {q} 
(tantque) si {I A C} A {1} alors 

{ 1} tantque C faire A ftantque {I A — iC} 


(I doit être un invariant de la boucle) 
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Exemple : règles d’affectation 

et conséquence 

• { p[ x / e] } x := e { p(x)} 

1. {n+1 > 0} n:= n+1 {n>0} (x=n, e = n+1) 

2. {k > 0} n :=0 {k > n} (x=n, e =0) 

3. {x = 4} x := x+1 {x=5} en écrivant: 

{(x+1 ) -1 =4} x := x+1 { x - 1 =4} => {x = 5} 

4. {x > 0} x := x + 1 {x > 1 } 

{x > 0} => {x + 1 > 1 } 

X := X + 1 
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Exemple : règles de condition 

5. [vrai} 

si x > 0 alors 

{vrai a x > 0} => {x > 0} =>{x+ 1 > 1 } 

X := X + 1 ; 

{x>l} 

sinon 

{vrai a x < 0} => {x < 0} => {-x > 0} 
x := - x; 

{x > 0} {x > 1 } 
fsi 

{X > 1} 
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Schéma de preuve d’un algorithme itératif 

Tant que 

□ Soit l’algorithme suivant: 

début 

{pré} 

init; 

{ 1 } 

tantque C faire 
{IaC} 

A 

{ 1 } 

ftantque 
{I a-,C}^ 

{post} 
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Exemple: calcul de la racine carrée 

par défaut 

■ Algorithme A: 
donnée : n : entier; 
résultat : r : entier; 
début 
r := 0; 

tantque n>(r+ l) 2 faire 
r := r + 1 ; 

ftantque 

fin 


- Spécification: 
précondition : {n > 0} 
postcondition : { (r 2 < n a n < (r + 1 ) 2 ) } 

■ invariant de la boucle: I = {r 2 < n} (évident, sinon on le déduit de la postcondition) 
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Algorithme A annoté et prouvé 


Algorithme A: 
donnée : n : entier; 
résultat : r : entier; 

Début 

{n > 0} => {n > 0*0} 
r := 0; 

{n > r*r} 

tantque n > (r + 1 ) 2 faire 

{n > r 2 A n > (r + 1 ) 2 } 
r := r + 1 ; 

{n > r 2 } 

ftantque 

{r 2 < n A n < (r + 1 ) 2 } 
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Exemple: 



■ Algorithme B: 
donnée : n : entier 
résultat : y : entier 
début 

X := n; 

Y := 1; 

tantque x > 1 faire 

y := y * x; 

X := X — 1 ; 

ftantque 

fin 


■ Spécification: 
précondition : {n > 0} 
postcondition : {y = n!} 

■ Invariant: I = { n! = y x! a x > 0} 


calcul de n! 
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Algorithme B annoté et prouvé 



Algorithme B: 
donnée: n : entier; 
résultat: y : entier; 

début 

{ n - 0 } 

X := n; 

{x > 0 A x = n} 

y := 1; 

{n! = y x! A x > 0} 

tantque x > 1 faire 

{n! = y x! A x > 0 A x > 1 } => {n! = (y x)(x-l )! A x > 1 } 
y := y * X; 

{n! = y (x-1 )! A x-1 > 0} 

X := X — 1 ; 

{n! = y x! A x > 0} => {n! = y x! A x > 0} 

ftantque 

{n! - y x! A x > 0 A x < 1 } => { y = n!} 

Exercice: Faites la preuve de l’algo. de n! en faisant une boucle qui va de 1 à n 
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